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ABSTRACT 

We report the detection of the (J,K) = a(4,3)--s(3,3) rotation-inversion transi- 
tion of ammonia at 124.6 /im toward the center of the Orion-KL region. The line is 
in emission and has a FWHM ^30 km s'^ The far-IR ammonia line emission prob- 
ably comes mainly from the "hot core", a compact region of warm, very dense gas pre- 
viously identified by the radio inversion lines of NH 3 . The a(4,3)— ^(3,3) line is very 
optically thick (t ~ 10 ^, and since it is seen in emission, radiative excitation of the 
(4,3) NH 3 level by far-IR emission from dust within the source can be ruled out. 
Radiative excitation via the 10 fim ro-vibrational transitions of NH 3 also seems 
unlikely. Hence, the (4,3) level is probably collisionally excited and the gas in the' hot 
core region is warmer than the dust. Since the far-IR line emission is highly trapped, 
densities of —10^ cm“^ are high enough to explain the observations. Shock heating 
by the mass outflow from IRc2 may account for the high gas temperatures in the hot 
core region^ 
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I. INTRODUCTION 

Interstellar ammonia was first discovered by its inversion transitions at 1.3 cm (Cheung et al. 
1968) and has become an important probe of the physical conditions in interstellar molecular clouds. 
Recent single dish observations of several of the 1.2 cm inversion transitions toward the center of the 
Orion-KL region have shown the existence of a component of "hot" ammonia gas in addition to the 
"spike" NH 3 emission from the quiescent molecular cloud (cf, Barrett et al. 1977; Wilson et al. 1979; 
Morris et al. 1980; Ziurys et al. 1981). High resolution mapping with the VLA has shown that the hot 
core NH 3 emission comes from a region of angular diameter about 10" (7 x 10*®cm) and temperature 
150 to 200 K. The opacities of many of the inversion lines are large (t ^ 10), and peak NH 3 column 
densities reach 5 x 10^® cm“^ (Genzel et al. 1982; Pauls et al. 1983; Palmer et al. 1983). The NH 3 hot 
core component may come from both streaming gas within and turbulent, swept up gas in interaction 
witliv outflow from the luminous infrared source IRc2 (Downes et al. 1981). Morris et al. (1980) have 
shown that the radio data are consistent with pure collisional excitation of the NH 3 levels if densities 
are 10^ cm”^, or with a combination of far-IR radiative excitation and collisional excitation at densities 
-“10^ cm~^ (assuming the NH 3 lines are optically thin). Intense emission at 3.5 mm from SO and SiO 
with characteristics similar to the NH 3 emission and from about the same region has been found by 
Plambeck et al. (1982) and Wright et al. (1983). Emission lines from vibrationally excited HC 3 N, 
CH 3 CN and torsionaily excited CH 3 OH have been detected in single dish observatiang of Orion-KL and 
probably also come from the hot core Orion-KL region (Clark et al. 1976; Loren et al. 1981; 
Goldsmith et al. 1982; Hollis et al. 1983). In the present paper, we report the first detection of a far-IR 
emission line of NH 3 . The new data are inconsistent with radiative excitation of the (4,3) NH 3 level, 
and support a collisional model. A hydrogen density of —10^ cm is all that is required, as the far-IR 
NH 3 lines are very optically thick. Collisional excitation at that density may also explain the emission 
of some of the other molecular lines from the hot core region. 
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II. OBSERVATIONS AND RESULTS ' 


The data were taken in February, 1982, with the 91.4 cm telescope on board the NASA Kuiper 
Airborne Observatory. The spectrometer wan a liquid helium cooled, tandem Fabry-Perot described by 
Storey, Watson and Townes (1980), with a photoconductive detector. The angular resolution was 44" 
(FWHM, or 55" equivalent disk), and the chopper throw was 4' at a chopping frequency of 29 Hz. The 
total system NEP (noise equivalent power) was 3 x 10~’‘* W Hz“*^. The J -"4— *3, a*-*s * NH 3 rota- 
tion inversion Lansitions at 125 /am in Orion-KL were observed at a resolving power X/AA 3400, 
resulting in a Lorentzian instrumental profile of FWHM 90 km s~^ Wavelength and velocity calibra- 
tion was provided by the NH 3 J '= 4 — * 3 lines in a gas cell, and an HDO line at 124.9547 /txm 
(McClatchey et al. 1973), with a precision of about ±10 km s“’. To calibrate the line intensities near 
125 m, we used the line to continuum ratios and assumed a flux density of 6.0 x lO'* Jy for the cen- 
tral 60" of the Orion-KL region at 125 /am (Werner et al. 1976). The J = 5— *4, a— *-s transition 
region at 100 /t,im was observed at a resolving power of X/ AX “ 4000, resulting in a velocity resolution 
(FWHM) of 75 km s“*. Absolute wavelength calibration was provided by the NH 3 lines in a gas cell, 
and a telluric H 2 ^^ line at 100.2601 /am. The observed 125 /am spectrum of Orion-KL is shown in 
Fig. 1, and the main results are as follows: 

I) The (J,K)“a(4,3)— *s(3,3) line (rest wavelength, 124.6474 ^m. Urban et al. 1981) is present in 
emission with an LSR velocity centroid of 0 ±10 km s"\ in agreement with the velocity centroids 
of the radio inversion lines (vlsr = 5 to 8 km s“‘). The observed line width (FWHM) is 
140±20 km s“\ which is significantly wider than the instrumental resolution (90±10 km s~*). 
Hence, the intrinsic velocity width of the line in Orion-KL has to be ^30 km s~^ if the line 
profile is Lorentzian. For a Gaussian shape the width would be substantially larger than this, and 
may be due to high velocity gas with too small a column density to have been detected in the 
inversion spectrum. The peak line flux (~ 5 x 10“*® W cm“^ within the pass-band of the spec- 
trometer) is ~ 6 % the intensity of the continuum at the resolution used. The effective line 

*s,a refer to the symmetry of the rotation inversion wave function with respect to reflection about the plane 
of the hydrogens. .\n a — > s transition is between the upper inversion level in the upper rotational state to 
the lower inversion level in the lower rotational slate, Each inversion level is further split by hyperfine struc- 
ture, which is not resolved in the present measurements. 
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brightness temperature is 25 to 30 K, assuming that the NH 3 sources fill the beam. 

2) The neighboring (J,K) = (4,2)— (3,2) line at 124.7957 /im (Urban et al. 1981) is not evident. At 
vlsr “ 5 km s“\ the 3 cr upper limit to the line intensity is 2 x 10“'® W cm"^ that is, about 3 
times weaker than the (4,3) line. There is a feature at Vlsr =“ +70 km s"' which, however, is 
not interpreted as significant, as its width is less than the instrumental resolution. 

3) The (4,1)— *(3,1) and (4,0)— *(3,0) transitions (rest wavelengths 124.8835 and 124.9125 /u,m, 
Urban et al. 1981) are close to the bottom of a telluric HDO absorption line. The offset to 
shorter wavelengths of the center of absorption in Fig. la from the HDO rest frequency, and what 
may be additional absorption at the wavelengths corresponding to vlsr ~5 km s“' is not incon- 
sistent with nonlinearities in the sweep due to distortion of the piezoelectric element in the 
Fabry-Perot drive. 

4) None of the observed J = 5— *4 transitions falling within observed wavelengths (K = 3, 100.1046 
/Ltm; K = 2, 100.2129 /x.m; K = 1, 100.2772 (xm. Urban et al. 1981) are present to a 3 o- limit 
greater than 5 x 10”'® W cm”^ (in absorption or emission). Only the (5,3)— *(4,3) transition is 
not confused by telluric absorption. The (5,4)— (4,4) line wavelength was not observed since this 
line was known to be almost coincident with a strong telluric H 2 O line. 

It may not be surprising that the (4,3)— (3,3) line is stronger than the other transitions, since it is 
the only observed transition whose upper level has K = J-1. Non- metastable levels with K<J-2 are 
expected to require substantially higher densities or a more intense radiation field to be populated 
equivalently to the K * J - 1 levels. 

III. DISCUSSION 

In the following, we discuss possible excitation mechanisms for the NH 3 gas and investigate which 
of the known gas components at the center of Orion-KL may account for the far-IR emission. Table 1 
lists the characteristics of the different components. In addition to the hot core and spike seen in the 
NH 3 inversion lines, there are also the "plateau" and "shocked gas" components in other molecular lines 
which probably come from gas within and in interaction with the "high velocity" outflow (e.g., Scoville 
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1981; Beckwith 1981). The spike feature has a linewidth of only a few km s \ and almost certainly 
cannot contribute signihcantly to the far-IR line (apart from lack of excitation, see below). 

1. Radiative Excitation 

The detection of an emission line immediately rules out that far-IR continuum radiation by dust 
mixed with the gas or by a far-IR source embedded in the line emitting region can alone account for 
the excitation of the (4,3) level. Far-infrared radiation by an internal source can only result in an 
absorption line or in a redistribution of the continuum radiation, but not in a net emission line. 
Resonant scattering by NH3 molecules behind or beside a continuum source could conceivably produce 
the emission line, but there is no observational evidence for a strong far-IR continuum source external 
to the center of the region (Wynn-Williams et al. 1983). Radiative excitation is also possible via the 
ro-vibrational bands of NH3, particularly the vj bands at 10 /am. Radiation from the luminous, compact 
source IRc2, for example, could be efficiently pumping the 10 /am transitions, since its spectrum peaks 
at about 8/am. However, detailed considerations, based on present knowledge of the structure of this 
region, rule out such excitation. To pump the far-IR transitions, the total number of 10 /am photons in 
a given ro-vibrational line should be greater than or equal to the observed number of photons in the 
far-IR line ((4,3) -> (3,3): N/-//? *= 1.5 ± 0.5 x 10‘*’s“*). The number of 10 ^m photons emitted by 
IRc2 per ro-vibrational transition is about = 5 x 10'*® A V/o (s“^), where Avio is the FWHM of a 
ro-vibrational line in units of 10 km s”*. This value is an upper limit for the number of photons avail- 
able for pumping the 10 /xm NH3 transitions in gas in the emitting region. For such a value, one must 
assume that the luminosity of IRc2 is 10^ L©, that is, equal to the luminosity of the whole Orion-KL 
region. Furthermore, the number of 10 ftm photons coming from IRc2 has to be corrected for emis- 
sion and absorption of dust along the path to the far-IR emitting region. Recent infrared observations 
suggest that there is a cavity of low dust density out to a radius of ~4 x 10^® cm from IRc2, which is 
surrounded by the dense, clumpy hot core region and a region which contains most of the quiescent 
and the high velocity gas (Werner et al. 1983; Wynn-Williams et al. 1983). The 10/am dust opacities 
through these regions are substantial (t ~3 to >10). Because of the presence of the large cavity 
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around IRc2, the dust grains in the hot core, spike, plateau and shocked regions absorb near infrared 
radiation from IRc2, but are too far from that source to reach temperatures so that they significantly 
emit at 10 /xm. The dust in the far-lR emitting region, therefore, cannot contribute significantly to 10 
/xm pumping. Since the far«infrared NH 3 lines are optically thick, the emission must come from the 
outer surface of the region, where the lOjixm radiation of IRc2 is attenuated by 1 to 4 orders of magni- 
tude. Therefore, the number of 10 ^m photons available in the far-lR emitting region, N//jc 2 x e“’‘, is 
probably much smaller than This probably rules out 10/xm radiative pumping. 

2 ) Collisionul Excitation: The Far-Infrared Emission Comes from the "Hot Core" 

We have used the values of temperature, hydrogen density and total mass within the field of view 
lusted in Table 1, to derive the maximum number of photons from (4,3) -> (3,3) transitions, for the 
different components of molecular gas in Orion-KL. For the hot core and spike, the column density of 
NHj can be directly estimated from the radio opacities. In other regions, the total number of NH3 
molecules is assumed to be approximately 10"^ that of Hj. Table 1 shows that only the hot core gas 
has sufficient density and numbers of NH3 molecules to produce the observed (4,3)— '(3,3) emission 
intensity. The spike region does not produce quite enough photons, and furthermore, its doppler velo- 
city width is much too small to provide the width found in the (4,3)— *(3,3) transition. Other regions 
would produce far too few collisional excitations. We propose, therefore, that the NH3 far-IR emission 
comes from the hot core region, and that the (4,3) level is colUsionally excited. As a consequence of 
the high opacity expected for the (4,3)— (3,3) transition from the hot core irpm ~600), the far-IR line 
also has to be broader than the 1.3 cm inversion lines, since amounts of gas at high velocity which are 
almost transparent at the inversion frequency can be observed in the rotational transition. This is con- 
sistent with the observed broadening of the 125 /xm line. If it is assumed that all of the 125 /xm con- 
tinuum comes from a source 10" in diameter, the radiation temperature would be about 200 K. An 
extrapolation of the 10 /xm opacity, with opacity proportional to a and a - 1 or 2, gives a 125 /xm 
opacity of 1 to 5, so this temperature should also be about equal to that of the dust. However, esti- 
mates based on the total luminosity and on the dust emission at 8 to 30 /xm give a dust temperature of 
^140 K for this 10" region. Some of the continuum radiation almost certainly comes from the 
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surrounding region and explains this discrepancy. The observed increase in emission at the 
(4,3)-*(3,3) transition of about 6% above the continuum intensity would be due to an effective tem- 
perature for the transition of 210 K if the entire emission comes from the 10" region. In any case, for 
the line to be in emission the rotational temperature must be higher than that of the dust; the above 
estimates show that the temperature difference is probably a few tens of degrees. 

The apparent absence of the (4,2)— (3,2) and (5,3)— (4,3) lines suggests that densities are not 
high enough to populate the (4,2) and (5,3) levels to the equivalent excitation temperature of the (4,3) 
transition. Microwave observations indicate that indeed the population of the (4,2) level is less than its 
thermal equilibrium value. The observed opacity of the (4,3) inversion line averaged over the source is 
consistent with thermal equilibrium at about 200 K (Genzel et al. 1982, Table 2). However, the 
observed opacity of the inversion line of (4,2) is between 0.5 and 0.9, about half that expected for 
LTE. While this is an indication that the (4,2) level is not in LTE, it is tlie population ratio between 
(4,2) and (3,2) states to which the far-IR is sensitive, rather than the population itself, since the gas is 
optically thick at 125 fim. For the (4,2) - (3,2) transition to show no emission at all, or perhaps an 
absorption of the dust continuum, its effective temperature must be reduced at least to 200 K and pos- 
sibly as low as 140 K. The latter, more stringent requirement implies that relative populations are 
(/i 4,2 “ « 3 , 2 ) / ('* 4,3 “ " 3 , 3 ) :^0-77. Thus, a 23% change in relative populations between upper and 
lower states due to a smaller ratio between the collisionai and radiative rates would explain the non- 
detection of emission from the (4,2)- (3,2) transition. This smaller ratio would have such an effect, of 
course, only if the excitation and radiative rates are comparable. For an optically thin medium, the 
radiative transition rates for the (4,3) and (4,2) states are proportional to (J* - through the rota- 
tional matrix elements. This gives a ratio for the two radiative rates of 0.58, with the (4,3)— (3,3) tran- 
sition more easily brought into temperature equilibrium by collisions than (4,2)— (3,2). Since the 
medium in fact has substantial optical depth, the smaller optical depth of the (4,2)— (3,2) transition 
further enhances the difference in excitation between it and the (4,3)— (3,3) transition. Table 2 gives 
parameters for several non-metastable NH3 levels. The values of rm itnd P(tfjr) are derived from 
the observed opacities of the 1.3 cm inversion lines. Values for collisionai cross sections between He 
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and NH 3 were taken from Green (1982), and were multipUed by 2.5 to account for the faster thermal 
speed of H 2 over He and for enhanced collisional cross sections of H 2 molecules which are in rotational 
states J > 0. The quantities and fix listed in Table 2 are the densities at which the popula- 
tion is half of that at LTE without and with trapping corrections, fi x no .9 is the density at which the 
population is 0.9 that at LTE, with trapping included. Table 2 shows that because of the large trapping 
corrections, densities of —10’ cm~^ are sufficient to account for the population of the (4,3) and (4,2) 
levels. Higher densities may be necessary to also account for the observed radio brightness tempera- 
tures of higher excitation lines, such as (7,6) and (10,9). Such high densities, however, may be incon- 
sistent with the weakness of the (4,2) radio and infrared lines. A possible solution might be a mixture 
of radiative and collisional excitation. At hydrogen densities of — lO’ cm~^ the NH 3 abundance in the 
hot core is — 10~^ to 10“^. A possible mechanism to heat the 1 to 10 Mqo( gas in the hot core region 
above the dust temperature may be the mechanical (shock) heating by the mass outflow from IRc2. 
The mechanical luminosity of this source is estimated to lie between 10 and 1000 L©, which is 
sufficient to account for the total line cooling from ammonia and other molecules. 
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Figure Captions 

Fig. 1. 125 /Am spectrum toward the core of the Orion-KL region, with a 44" beam FWHM. The 
spectral resolution is a Lorentzian of FWHM 90±10 km s“l 

Lower (a): Observed spectrum, with arrows marking the positions of the a(4,K)~*s(3,K) rotational 
inversion transitions of NH 3 toward Orion (vlsr " + 5 km 8 “’). Also marked is a telluric HDO line 
on the left side of the spectrum. The K — 0 and K — 1 NH 3 lines are close to the bottom of this 
absorption feature. There is a second HDO absorption line at 124.3008 /tm which causes the downward 
slope on the right side of the spectrum. 

Upper (b): The same spectrum, but smoothed and with Lorentzians fitted to the HDO lines and then 
divided out. The positions of the K * 2 and K — 3 NH 3 lines are indicated, together with the 
appropriate velocity scales. A Lorentzian fit to the (4,3)—*(3,3) transition (thin line) gives an LSR 
velocity centroid of 0±10 km s~\ and a FWHM of 140 ±20 km s“', that is, significantly wider than the 
resolution. The K — 2 line is at least three times weaker than the K — 3 line. The positive bump at 
LSR + 70 km s~* is probably narrower than our resolution, and therefore is not interpreted as real. 
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t Tlie niaximuin FIR pholons/sec whicli can be ciiiiUcd or absorbed is taken to be 
Nn.i, iifi a- V exfil— h p I kT) where Nnh^ is the number of NHj molecules in the J — 3, K “ 3 state if 
xNHj- I0~’; is assumed to be 10"*®cm^s"*, and ftp corresponds to the rotational energy difference. 
This is an upper limit to the number of photons emitted if the line has substantial opacity. 
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Collisional Excitation of the NH^ Levels 
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The opacities of the 1.3 cm inversion lines are estimated from the single dish observations by Morris et al. 
1980, Zuckerman et al. 1981, Ziurys et al. 1981, and represent an average value across the source. For the 
(4,2), (7,6) and (10,9) lines, the opacities are obtained by comparing to the line temperature of the optically 
thick (4,3) line. 

The departure coefficients b •» nupptr! '^upptr (LTE) are estimated from the observed opacities and the opa- 
cities extrapolated from r(3,3) * 20 at T = 200 K. 

1 — 

The escape probability for far-infrared line radiation: /8(r) - — — . The value used is a compromise 
between an expanding source and Gaussian turbulence (de Jong et al. 1975, 1980). 

is the density where 6 - without trapping of FIR radiation; 
(b = nupper! ^upperUTE) * 0 + ^ui C„,)"', where is the Einstein coefficient and C^i the colli- 

sional rate between upper and lower level). 

I3xn^^ is the density where b — 'k if trapping is taken into account, 
is the density where b - 0.9 if trapping is taken into account. 








